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Abstract 

Isospin symmetry breaking in the quark condensates, (dd) ^ {uu), is a fundamental 
parameter in both the QCD sum rule studies and the chiral perturbation theory. In 
this article, we apply the QCD sum rule method to treat the hyperon mass difference 
i(m£- + 777.£o) — m^+ in order to obtain a reasonable estimate on such isospin symmetry 
breaking parameter 7 defined as ({dd) — {uu))/{uu). Note that the electromagnetic 
contributions to the particular mass difference are expected to cancel almost completely. 
Using the light-quark mass difference — m u = 4 ± 1 MeV and the experimental data 
^( m s- + "^s°) — "i£+ = 5.62 ± 0.13 MeV, we obtain 7 = —0.011 ± 0.001, a value which 
is slightly larger than the commonly adopted value and is of smaller uncertainty. 



PACS number(s): 11.30.Hv; 12.38.Lg; 11.55Hx; 14.20.-c 



1 



Although quantum chromodynamics (QCD) has been accepted as a theory for describing 
strong interactions among quarks and gluons, it is highly nonperturbative at the confinement 
scale. In particular, it is expected that the QCD ground state, or the QCD vacuum, has 
nonzero quark condensate < qq >, gluon condensate < G 2 >, and many other condensates of 
higher order. Making use of the method QCD sum rules]]]], it has become possible to study 
hadron physics in terms of the basic parameters of QCD, such as current quark masses and non- 
perturbative condensates. In a systematic study of isospin symmetry breaking effects, both 
the quark mass difference, m d — m u , and the difference in condensates, 7 = ((dd) / (uu)) — 1 are 
the most fundamental parameters, while the quark mass difference — m u can be estimated 
with reasonable confidence as — m u ^3 — 5 MeVp], [!|, 7 is poorly known. Theoretical 
estimates on 7 have been obtained in the chiral perturbation theory [Q, in the Nambu-Jona- 
Lasinio model[]5], and in QCD sum rules[0, [| ^ [10], |TTJ; the value has a fairly large range, 
7 = — (2 ~ 10) x 10 -3 . In the QCD sum rule approach, the error in estimating the value of 7 



comes from the poorly-known electromagnetic contributions ]8|] and instanton effects [|TT . 

In this paper, we wish to use the QCD sum rule method to obtain the isospin symmetry 
breaking parameter 7 from the hyperon mass difference ^(m^- + mjo) — m-£+, a combination 
that the electromagnetic contributions are expected to cancel almost exactly. The cancelation 



may be understood as follows. In the simple (constituent) quark model []12|, [13], y, |3j, the 
electromagnetic contributions can be separated into a contribution from photon exchange 
between two of the three quarks and a contribution from the quark self-energy type: 

= a(l/r)(Q 1 Q 2 + Q 2 Q 3 + Q3Q1) + C(Q? + Qj + Q 2 3 ), (1) 

where Qi, Q 2 , and Q 3 are the quark charges and (1/r) is the (effective) average inverse distance 
of two quarks in a baryon. In the isospin symmetric limit, the radii of S + , E~, and S° are 
identical so that 

(l/r) s+ » (l/r) E - « (l/r) s o, 

C S+ ^C S "^C S °. (2) 



Here the equality would take hold should isospin symmetry be exact. We therefore obtain 

-Os- +ml ) -ml + = 0. (3) 

We note that similar arguments may be constructed for other models of hadrons. Using 
the quark model estimates and the experimental information from electron scattering on the 
protons and neutrons (extrapolated to the other members of the baryon octet), Gasser and 
Leutwyler0 have reported |(?7t£- + m s°) — m E+ = — 0.32±0.26 MeV, which is consistent with 
the above observation. Accordingly, the experimentally observed value, i(ms-+mso)-ms+ = 
5.62 ±0.13 MeV, may serve to determine m^ — m u and 7 in an unambiguous manner. As some 
authors have indicated that instanton contributions are absent in the chirally even baryon 
sum rules|ll[] but may be of significance in the chirally odd one, we choose to employ solely 
the chirally even sum rules. 

We turn our attention to the derivation of QCD sum rules. Following the standard proce- 
dure, we consider the two-point Green's function 

n(p) = ifd 4 x e ipx (0|T(r/(a;)?7(0))|0), (4) 

where the composite field operators are specified by 

= e ahc ^=[(u aT C^d b ) l5 ^s c + (d aT C^u b )^s% 

Vs - = e abc (d aT C^d b ) l5 Ys c , (5) 

with a, b, and c color indices and C=—C T the charge conjugate operator. At the hadronic 
level, we write the Green's function via dispersion relation as 

„. . 1 r°° Im U(p') , , 9 

7T JO P U — P Z — IS 

where 

Im n(p') = nJ2$(p' 2 - M$)(0\v(0)\X)(X\m\0), (7) 
x 



and 

(Ofo(0)|E) = A E «(p) > (8) 

with the Dirac spinor normalized as v(p)v(p)=2m-£. At the quark level, we write, using 
the S + case as our example, 

(mv{x)nm\V)=^ abc t a ' b ' c 'Tr[Sf (9) 

Here we have used the definition, 

zS ab (x) = (0\T[q a (x)q b (0)]\0), 

as in Ref. §. 

has two invariant structures: 

u( P ) = u 1 ( P 2 ) + i^( P 2 ) a (10) 

where 111 and 11^ are two invariant functions. We name 111 as the chirally odd sum rule 
since chiral-odd operators contribute dominantly, and ITy as the chirally even sum rule. The 
derivation of the relevant QCD sum rules is similar to what we have reported in our earlier 
papers [§, I|Jj. Here we record only the final results on the H^ (chirally even) £ mass sum 



rules: 



32 6 24M 2 



m Q a M 



£-4/9^£+ _ '"'s^siii'O ^-26/27 _ '"'"""'"O ^-26/27 



I ' " !l ~ 24 ~~ 12 



,o2 -(M^/Af 2 ) = ^_ r-4/9 77S , ^ r-4/9pS° , Q tt Q rf ^4/9 _ a « a d m Q r-2/27 

Ps ° 8 2 32 6 ~ 24M 2 

[m s a s - (a u - a d ) (m d - m u )]M 2 _ 4/9 s o m s a s m' 2 



_£-4/9^£° _ "^^'"O ^-26/27 



and 



4 24 

[flrf(3m« - mj + an(3m d - m n )]mg r _ 26/27 , . 

48 ' 1 j 



2 (M 2_ /M2) = 4/9 E - b_AP 4/9 x- 4^/9 _ « £-2/27 

Ms 8 2 32 6 24M 2 



2 



— — — L _4//9 _E S_ — msasm ° j^-26/27 _ m rf Q rf m Q ^-26/27 



4 " 24 12 

where we have adopted the definitions: 



a 9 = -(2vr) 2 (gg), 



? mg = (2n) 2 (g s qa ■ Gq), q = u or d, 
a s = -(27r) (ss), 



a s md = (2rr) (g s s<7 • Gs). 
$ - (2vr) 4 ^, 



= 1 - e-*c»o(l + x (m) + . . . + ^xf m) ), with x (m) = Wf m) / M 
_ _ ln(M 2 /A 



r2 / A 2\ 



■ 



(14) 



ln(u 2 /A 2 ) 

In Eqs. (11-13), M 2 is the Borel mass squared, as arising from the Borel transform of the sum 
rules, 

B[U m (p 2 )}^ lim ^(-P 2 )" +1 (^) n n lW (p 2 ). (15) 

— p 2 — >oo ' " 

~ p a fixed 
nM 2 J 

In what follows, we take the value (ss) =0.8 x((uu) + (dd))/2. We also use the relation 

(g s scx ■ Gs) / (g s ucr ■ Gu) ~ (ss) / (mm) , 

which indicates mj 2 ~ m 2 ,, a relation which has been examined in some detail in octet baryons 
J/J. The input parameters in the numerical analysis are taken to be a q =(a u + ad)/2=0.545 
GeV 3 , a s =0.8 a q , 6=0.474 GeV 4 , m§=0.6-0.8 GeV 2 , m q =(m u + m d )/2= 7 MeV, m d - m u =3-5 
MeV, and m s =150-200 MeV. These values correspond to the QCD scale parameters A= 0.1 
GeV, and the normalization point «= 0.5 GeV. Applying the operator M 4 d£n/dM 2 to both 
sides of Eqs. (11-13) (to get rid of the dependence on the parameter /3f,), we deduce the QCD 
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sum rule on the mass difference |(j71e- + m E o) — m^+. We fit the two sides of this sum rule 
by minimizing 

S(M 2 ) = ]T \LHS - RHS\ 

M 2 =1.0+0.005i 
0<i<60 

within the Borel region 1.0 GeV 2 < M 2 < 1.3 GeV 2 , where the LHS (the left-hand side of the 
sum rule) is the experimental value, while the RHS (the right-hand side of the sum rule) is 
the operator-product-expansion results at the quark-gluon level. In the numerical analysis, we 
adopt W s+ =3.47 GeV 2 flj and the experimental data |(m s - +m s o) -m s+ =5.62±0.13 MeV. 
Our results are (W 2 + W/ 2 „)/2=3.483±0.002 GeV 2 , W 2 ^ - 1U 2 =0.0030±0.0004 GeV 2 , and 
7 = —0.011 ± 0.001. Theoretical uncertainties in our determination of the isospin symmetry 
breaking parameter 7 come mainly from the uncertainty in the value of — m„ and the 
experimental error in |(mg- + m^o) — m-^+. 

In Fig. 1, we plot the LHS and RHS as a function of the Borel mass squared, M 2 . The 
solid curve is the LHS assuming the value n( m s- + m E°) — m s +=5.62 MeV. The dotted curve 
is the RHS with the optimized values, 7 = -0.0114, W| =3.4816 GeV 2 , and W 2 . =3.4846 
GeV 2 , where we have used, as the input, m = 0.8 GeV 2 , m s = 175 MeV, md=8.9 MeV, and 
m u =5.1 MeV. 

Following the same procedure as outlined above and using the present value for 7, we 
have repeated the analysis of the sum rule for the neutron-proton mass difference within 
the Borel region 0.9 GeV 2 < M 2 < 1.1 GeV 2 @, |TJ. The result is illustrated in Fig. 2, 



in which the solid curve is the LHS assuming the value m^ D — rnji CD =2.693 MeV and 
the dotted curve is the RHS with the optimized values, 7 = -0.0102, 14^=2.2506 GeV 2 
(neutron), Wp=2.25 GeV 2 (proton), rrid=S.9 MeV and m„=5.1 MeV. Such value implies that 
the electromagnetic contribution to the neutron-proton mass difference m? N — m? p is about 
— 1.3 ~ —1.8 MeV, a value in reasonable agreement with the naive quark model result P|. 
Note that our result on 7 is a little larger than that obtained by Forkel and Nielsen []TT|, who 
obtained 7 = —(0.8 — 1.0) x 10 -2 making use of rn' N — m? P =— 0.76±0.3 MeV[Q (corresponding 
to i.e., m% CD -m Q P CD = 2.05±0.3 MeV) in the Borel region 0.8 GeV 2 < M 2 <1.4 GeV 2 . Using 

6 



their value on mjf — m? P and adopting our approach together with the same parameters as in 
the present paper, we obtain 7 = —(0.7 — 0.9) x 10™ 2 , a result consistent with that in Ref.[T 1 



and also with that in our earlier paper ||. Nevertheless, the choice rn^ — irtp = —0.76 ± 0.3 
MeV can only be marginally consistent with the result obtained through the analysis of the 
sum rule for |(m s - + m s o) — m s +, which gives a better determination of 7 and leads to 
— nip is about —1.3 ~ — 1.8MeV. 

In summary, we have presented a detailed analysis of the isospin symmetry breaking pa- 
rameter in the quark condensate, 7, from the chirally even QCD sum rules on the mass 
difference + m s o) — m^+, for which the instanton contributions are absent and the 

electromagnetic contributions are expected to cancel almost exactly. Our numerical analysis 
indicates that a value of 7 (= —0.011 ± 0.001) (slightly larger in magnitude) is preferred over 
the existing result, 7 = —(2 ~ 10) x 10~ 3 . Using the present value of 7 and the observed 
neutron-proton mass difference, we have also deduced the electromagnetic contribution to the 
neutron-proton mass difference as in the range —1.3 ~ —1.8 MeV, in reasonable agreement 
with the quark model result. 
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Figure 1: |(m s - + m s o) — m s + as a function of the Borel mass squared M 2 . The solid curve 
is the LHS, the left-hand side of the sum rule, assuming |(m s - + m s o) — m s +=5.62 MeV. 
The dotted curve is the RHS, obtained by making use of m = 0.8 GeV 2 , m s = 175, m d =8.9 
MeV, and m u =5.1 MeV as the input. 
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Figure 2: mjf — ml as a function of the Borel mass squared M 2 . The solid curve is 
the LHS assuming m% CD — rrip CD =2.693 MeV and the dotted curve is the RHS with the 
optimized values, 7 = -0.0102 and W#=2.2506 GeV 2. Here we have used W^=2.25 GeV 2 , 
m^=8.9 MeV and m u =5.1 MeV as the input. 
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